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ABSTRACT 



The Ising spin glass model in a transverse field has a zero temperature phase 
transition driven solely by quantum fluctuations. This quantum phase transition 
^ . occuring at a critical transverse field strength has attracted much attention re- 

cently. We report the progress that has been made via Monte Carlo simulations 
of the finite dimensional, short range model. 

The most frequently discussed spin glass models are classical systems for which 
quantum fluctuations can be neg lected@. In most cases this is correct, namely as long 
as T c > 0. The reason for that is that critical fluctuations at the transition occur at 
a frequency u c (= r _1 , where r is the relaxation time) with Tiuj c <C ksT since u c — > 
as T — > T c due to critical slowing down. Hence any finite temperature will destroy 
quantum coherence and the system will behave classically. Very recently however, 
spin glasses began to enter the quantum regime!. 

The interesting theoretical question is: What are the effects of quantum mechan- 
ics on the physics of strongly disordered systems at zero temperature, where no heat 
^ ■ bath is present and hopping over energy barriers is replaced by tunneling through 

them quantum-mechanically? The renewed interest in spin glasses in the quantum 
regime was kindled by a series of recent experimentsi on the dipolar Ising magnet 
Li a .Hoi_ a; YF4, where T c was driven down to zero by the application of a transverse 
magnetic field T (see the phase diagram depicted in figure [I]). Therefore it became 
possible to study the zero temperature phase transition occuring for critical transverse 
field strength T c by simply tuning the external field. This quantum phase transition 
lies within a different universality from the usually studied (classical) spin glass tran- 
sition at finite temperatures and it turns out that their properties differ significantly. 

The above mentioned experiments can be described by the model Hamiltonian of 
an Ising spin glass in a transverse field^l 

h = -Y.J^ -r£< > (!) 

(ij) i % 

where the <7j are Pauli spin matrices, T is the strength of the transverse field and h is 
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Fig. 1. Typical phase digram of a three- (left) and two- (right) dimensional quantum Ising spin 
glass in a transverse field. 

a longitudinal magnetic field used to define magnetic susceptibilities but usually set 
to zero. Obviously, for T = the quantum-mechanical Hamiltonian ([I]) is diagonal in 
the z-representaion of the spin operators, which in this case can simply be replaced 
by their eigenvalues ±1 (after rescaling the couplings) giving exactly the classical 
EA-model in d dimensions: H = —J2(ij) JijSiSj, where now Si — ±1 are Ising spin 
variables. In this way the transverse field introduces quantum mechanics into the 
spin glass problem and the value of T tunes the strength of the quantum fluctuations. 
At zero temperature and T = the system described by (|1|) will be in its uniquely 
determined ground state, which is identical to the classical ground state of the EA- 
spin glass model. In this case one has (erf) = ±1 for all sites i and therefore qEA = 
K°f ) 2 ]av = 1 ■ where (• ■ ■) means the quantum-mechanical expectation value. 

If we switch on the transverse field (r > 0) the Hamiltonian (|1|) is not diagonal 
in the ^-representation any more and its ground state will be a superposition of 
the classical ground state plus various excited states, which describes the quantum- 
mechanical tunneling at zero temperature between the local energy minima of the 
classical Hamiltonian. Furthermore |(of)| < 1 since the transverse field tries to align 
the spins in the x-direction and therefore qEA — [(°f) 2 ]av < 1- Increasing V diminishes 
the EA-order parameter qEA and for some critical value T c it will be zero: qEA = 
for T > T c . This is the zero temperature phase transition we are interested in and 
obviously we cannot expect that its critical properties have anything in common with 
the finite temperature classical spin glass transition. 

In order to describe this zero-temperature transition one introduces a quantity 
measuring the distance from the critical transverse field strength (at T = 0) 5 = 
(r — r c )/T c . If one assumes a conventional second order phase transition one has 
Qea ~ l^l' 3 and spatial correlations decay on a characteristic length scale that diverges 
at the critical point as usual: £ ~ \S\ ~ v and these exponents defined so far would be 



sufficient to describe the static critical behavior of a classical spin glass transition 
(provided hyperscaling holds). However, at a zero temperature transition driven 
solely by quantum fluctuations static and dynamic quantities are linked in such way 
that the introduction of a characteristic time-scale (or inverse frequency) is necessary: 
£ T ~ £ z ~ |<5|~ zz/ , where z is the dynamical exponent. This will become more evident 
below, when we consider an equivalent classical model. 

Much work in the past has been devoted to the infinite range modefl and the 
phase-diagram in the T-T-plane looks similar to the one shown in figure [l], i.e. for low 
enough temperature T and field T one finds a transition line separating a paramag- 
netic phase from a spin glass phase. Recent work! focused on the zero-temperature 
critical behavior and calculated the critical exponents 7 = 1/2 (with multiplicative 
logarithmic corrections), (3 = 1 and zv = 1/2. 

In two and three dimensions, which are most relevant for the above mentioned 
experiments, no analytical results are known. For this reason extensive Monte Carlo 
simulations have been performed recently in two dimensions by Rieger and Young! 
and in three dimensions by Guo, Bhatt and HusJ. Usually the investigation of 
quantum systems via Monte Carlo methods are hampered by various deficiencies, 
the sign problem being the most notorious one in this respect (see§ for a review). 
In studying the Ising spin glass in a transverse field however, one can exploit the 
fact that it can be mapped exactly onto a classical Ising model described by a real 
Hamiltonian. Using the Suzuki- Trotter formula! one can easily show that the ground 
state energy of the <i-dimensional quantum mechanical model (|1]) is equal to the free 
energy of a (d+1) -dimensional classical model, where the extra dimension corresponds 
to imaginary time, i.e. 

E(T = 0) , T ,„ 1 1 m o 

i— 5 — '- = lim — 7 Tr e~ pH = -rTr e" 5 2 

L d T^oL d At L t L v ; 

where the imaginary time direction has been divided into L T time slices of width At 
(AtL t = (3), and the effective classical action, S, is given by 

5 = - E E KtjSiWSjiT) -EE mwSiiT + 1) - E E ^(r) , (3) 

T (ij) t i r i 

where the Si{r) = ±1 are classical Ising spins, the indices i and j run over the sites of 
the original <i-dimensional lattice and r = 1, 2, . . . , L T denotes a time slice. Moreover 
it is Kij = ArJy, H = Ar/i and exp(— 2K) = tanh(Arr). One has the same random 
interactions in each time slice. In order to fulfill the second equality in @ precisely, 
one has to perform the limit At 0, which implies — > and K — > oo. However, 
the universal properties of the phase transition are expected to be independent of At 
so we take At = 1 and set the standard deviation of the to equal K. Thus K, 
which physically sets the relative strength of the transverse field and exchange terms 
in ([[]), is like an inverse "temperature" for the effective classical model in (|3|). 



One sees that the (d+l)-dimensional classical model (|3]) should order at low "tem- 
perature" (or coupling constant K) like a spin glass in the d spatial dimensions and 
ferromagnetically in the imaginary time direction. From this one concludes the exis- 
tence of two different diverging length-scales in the classical model ([3]): one for the 
spatial (spin glass)-correlations, which is £, and one for imaginary time (ferromag- 
netic) correlations, which is £ T . Thus in the representation (||D the link between statics 
and dynamics in the original quantum model ([!]) becomes most obvious. Correspond- 
ingly, to analyze the critical properties of the extremely anisotropic classical model 
(HD one has to take into account these two length scales via anisotropic finite size 
scaling^. 

Monte Carlo simulations of the classical model (§) are straightforward but it 
turns out that sample-to-sample fluctuations are significant, so one has to do an 
extensive disorder averagei'l. However, the finite size scaling analysis is complicated 
by the existence of two diverging length scales £ and £ r and one has to deal with 
two independent scaling variables: as usual L/£ and in addition the shape (or aspect 
ratio) L T /L Z of the system0. Thus, with the usual definition of a spin glass overlap 
Q = L~ d L^ 1 J2i, T Si(. T )Si(. T ) ( a an d b are two replicas of the same system) for the 
classical system, the dimensionless combination of moments of the order-parameter 
g av obeys 

g av {K, L, L T ) = 0.5[3 - (Q 4 )/(Q 2 )\ V ~ g w {8L l ' v , L T /L Z ) . (4) 

Here (• • •) means the thermal average and [• • -] av means the disorder average. In 
isotropic systems one has z — 1, which makes the aspect ratio constant to one for the 
choice L = L T and in order to determine the critical coupling K c one exploits the fact 
that g av (K, L, L) becomes independent of L for K = K c . In the present case of a very 
anisotropic system z is not known a priori and one has to vary three different system 
parameters to obtain an estimate for K c and z (and other exponents). The following 
methodic enhances the efficiency of such a search in a three-parameter space and 
also produces reliable estimates for the quantities of interest: In the limit L T ^> L z 
the classical (d + l)-dimensional classical system is quasi-one- dimensional, and in 
the limit L T <C L z the system is quasi-<i-dimensional and well above its transition 
"temperature" in d dimensions (which is even zero for d = 2). Therfore one has 
g aY {x,y) — ► both for y — >• and for y — > oo. Hence, for fixed x, gg, v (x,y) must 
have a maximum for some value y = y m & x {x)- The value of this maximum decreases 
with increasing L in the disordered phase K < K c (where 5 = (K c /K — 1) > 0) 
and increases with increasing L in the ordered phase. This criterion can be used to 
estimate the critical coupling, as exemplified in figure (0). If one plots g av (K c , L, L T ) 
versus L T /L Z with the correct choice for the dynamical exponent z, one obtains a 
data-collaps for all system sizes L. Finally one uses systems with fixed aspect ratio 
L T /L Z to determine critical exponents via the usual one-parameter finite size scaling. 
Various scaling predictions can be made if one supposes a conventional second 
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Fig. 2. The dimensionless cumulant g av for the (2+l)-dimensional classical model (0) for three 
different values of the coupling constant versus the system size L T in the imaginary time direction. 
The system size in the space direction is L = 4 (o), 6 (+), 8 (□), 12 (x) and 16 (A). Since the 
maximum of g av (L T ) is roughly independent of L at K^ 1 w 3.30 one concludes that the latter value 
is the critical coupling constant. 

order phase transition to occur at some critical "temperature" K c for the classical 
model. Let us assume that we are at the critical point K = K c and the aspect 
ratio stays constant L T oc L z . Then the order parameter scales with system size as 
Qea = [\(Q)\]av ~ L~PI V and the corresponding susceptibility xsg = L T L d [(Q) 2 ] av ~ 
IS i ' v with 'j/u = d + z — 2/3 /v. The equal time correlation function in the infinite 
system decays like C(r) = [(5 , j(r)S'j +r .(r)) 2 ] av ~ r -{<i+z-2+ ri ) ^ f rom w hich one obtains 
by integrating over r the scaling relation 2/3 /v = d + z — 2 + 7/. The (imaginary)- 
time-dependent autocorrelation function of the infinite system decays like G(t) = 
[{Si(r)Si(t + r))] av ~ j : ~{d+z-2+r))/2z an( j integrating over t one gets the uniform 
susceptibility Xf(w = 0) = d[(of)] av /<9/i ~ J_risl v with 7^ = z — (3/v. Thus Xf wm 
diverge if z > /3/v. 

The nonlinear susceptibility Xni = d 3 [{a z )] SiV /dh 3 ~ lr"tnil v [ s m ost relevant for 
the experiment si. This quantity, also given by = [3(M 2 ) 2 — (Af 4 )] av /L d , with 
M = Yji, T Si(r), can be shownB to scale in the same way as L d L 3 [(Q 2 )] av , which 
yields the scaling relation 7^ = v{2 — 77 + 2z). The nonlinear susceptibility of the 
infinite system obeys Xni(T,8) ~ 5~ lnl Xni(T / 5 zu ^j with Xnz(^) ~~ > const, for a; — ► 
and Xni(x) — > x~" /nl ^ zu for x — > 00. By setting T = T c one gets Xnz ~ T^ ln ^ zu and 
similarily ~ r £~* i fl xv . 

In table [I] we list the results obtained so far in various dimensions. The symbol 
div. in the second column means that first and higher derivatives of the magnetiza- 
tion diverge already in the disordered phase. The word finite means that in three 
dimensions and for d larger than the upper critical dimension the uniform suscepti- 
bility does not diverge at the critical point. The results in the second column@ were 
obtained within a renormalization group calculation, those in columns 3 to 5 with 
Monte Carlo simulationsEl'ifl, column 6 shows the result of a Migdal-Kadanoff RNG 



calculation and the last column depicts analytical results from mean-field theorytl. 
Note that in the case d = 1 we refer to the non-frustrated random transverse Ising 
chain, where the order parameter is the magentization instead of the spin glass overlap 
and the exponents listed refer to the critical behavior of quantities that are defined 
accordingly. 



Table 1. Summary of critical exponents characterizing the quantum phase transition of the transverse 
field Ising spin glass in finite dimensions 
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These results imply in three dimensions a strong divergence of the nonlinear sus- 
ceptibility when approaching T = 0: Xni ~ T~ 2 - 7 (since j n i/isz ~ 2.7), with a power 
that is very close to the one for classical Ising spin glasses (~ 2.9, reffl). This is in 
striking contradiction to the observation made in the experiments mentioned in the 
introduction! and in order to clarify this discrepancy further work is necessary. 
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